Simulation of high-spin Heisenberg models in coupled cavities 
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We propose a scheme to realize the Heisenberg model of any spin in an arbitrary array of coupled 
cavities. Our scheme is based on a fixed number of atoms confined in each cavity and collectively 
applied constant laser fields, and is in a regime where both atomic and cavity excitations are sup- 
pressed. It is shown that as well as optically controlling the effective spin Hamiltonian, it is also 
possible to engineer the magnitude of the spin. Our scheme would open up an unprecedented way 
to simulate otherwise intractable high-spin problems in many-body physics. 



The Heisenberg spin model has played a crucial role as 
a basic model accounting for the magnetic and thermo- 
dynamic natures of many-body systems. Despite exten- 
sive investigations, however, many aspects of the model 
are still largely unexplored both analytically and numer- 
ically, especially for the cases of higher spins. The main 
difficulty in the numerical treatment originates from the 
fact that the Hilbert-space dimension blows up exponen- 
tially as the number of spins increases. As Feynman first 
noted [l|, this difficulty would be overcome in terms of 
quantum simulation based on precisely controlled quan- 
tum systems. Realization of quantum simulation, ex- 
pected in a near future, will mark a milestone towards 
the realization of sophisticated quantum computation. 

In the context of quantum information processing, a 
qubit is identical to an s — ^ spin, and in a few im- 
plementations, such as the arrays of Josephson junc- 
tions [2] or quantum dots 3|, the spin-chain Hamilto- 
nian naturally emerges from the spin-like coupling be- 
tween qubits, albeit with limited control of the coupling 
constants. On the other hand, in optical lattices, per- 
turbative evolution with respect to the Mott-insulator 
state can be described by an effective spin-chain Hamil- 
tonian 0, [j| . This approach has its own merit in that the 
spin-coupling constants can be optically controlled to a 
great extent. An alternative approach, recently under ac- 
tive investigation, is to use the array of coupled cavities, 
which are ideally suited to addressing individual spins 
[1, 0, 1 !, M M El El M- In th is approach, a spin 
is represented by either polaritons or hyperfine ground 
levels. The former, proposed in Refs. [H and allows 
a stronger spin-spin coupling than the latter, but lacks 
the optical control of the coupling. On the other hand, 
the latter, proposed in Ref. [9j, retains the optical con- 
trollability, but relies on rapid switching of optical pulses 
and the consequent Trotter expansion, which unavoid- 
ably involves additional errors and makes error-free im- 
plementation more difficult. More importantly, the ques- 
tion of simulating chains of higher spins, which may have 
a completely different phase diagram, remains open. In 
some sense, these are more important to simulate because 
unlike spin-i chains, they do not have exact analytical 
solutions for a wide range of parameters including the 



XXX case, except for special kinds of models [la ]. Addi- 
tionally, going to higher spins should make perturbative 
spin-wave theory more accurate, whose predictions can 
be tested. 

In this paper, we propose a scheme to realize the 
anisotropic (XXZ) or isotropic (XXX) Heisenberg spin 
model of any spin in an arbitrary array of coupled cavi- 
ties. Our scheme is experimentally feasible in that simply 
applying a small number of constant laser fields suffices 
for our purpose. If the number of lasers is increased, 
the individual constants of the spin Hamiltonian are con- 
trolled more flexibly. Most of all, a strong advantage of 
our scheme is that the magnitude of the spin itself can 
be engineered arbitrarily. This advantage contrasts with 
all the earlier schemes mentioned above including those 
for optical lattices, in which the spin is fixed in nature 
mostly as s = | (s = 1 in Ref. [||). s > \ spin chains 
exhibit fascinating physics that s — ^ spin chains can not 
have. A well-known example is Haldane's conjecture that 
antifcrromagnetic Heisenberg integral-spin chains have a 
unique disordered ground state with a finite excitation 
gap, whereas half-integral-spin chains are gapless [l6l.[l7|. 
Our scheme could be used to prepare a ground state, for 
example, through an adiabatic evolution, and measure its 
excitation gap and spin correlation functions [|[ . The ad- 
vantage of our scheme is, however, more apparent when 
we consider higher-spin problems intractable with any 
previous method. One of the intriguing examples is the 
quantum spin dynamics of the ferric wheels such as Feg 
and Feio composed of s = | spins [l8|. These systems 
would be simulable with a relatively small number (6 or 
10) of cavities. Spin chains also play an important role as 
a quantum channel for short-distance quantum commu- 
nication [19( . The property of an s — I antifcrromagnetic 
spin chain as a quantum channel strongly depends on its 
phase [2(J. In some phases, it provides an efficient chan- 
nel, outperforming that of a ferromagnetic chain. It has 
been also shown that a ground state with an excitation 
gap, as is the case for the spin-1 chain, can serve as a 
more efficient quantum channel plj . The ground state of 
the antifcrromagnetic spin- 1 chain also establishes the so- 
called localizable entanglement between two ends, which 
can be extracted by measuring every intermediate spins 
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FIG. 1: Involved atomic levels and transitions. Both transi- 
tions | a) <-> |e) and \b) <-» |e) are coupled to the same cavity 
mode with coupling rates g\ and gi and with detunings Ai 
and A2, respectively. Two laser fields with Rabi frequency ill 
and f^2 are also applied with detunings Ai and A2, respec- 
tively, and the transition between ground levels \a) and |6) is 
driven with Rabi frequency cu/2 by Raman lasers. 7 denotes 
the atomic spontaneous decay rate. 



[22| and then used for quantum communication. It is 
hard to demonstrate these schemes in optical lattices ow- 
ing to the difficulty of addressing individual sites [23| . 
Although the idea of communicating using spin chains 
is ultimately meant for solid-state applications, our sys- 
tem can serve as a preliminary test for comparing and 
contrasting the performance of various spin-s chains. 

We use two ground levels of a three-level atom to rep- 
resent ans=| spin (in a rotated basis, as will be seen 
later). We start by recalling that in terms of two states 
||) and |t) of one atom, the s — ^ spin is described in 
terms of operators s z = ±(|T) (T| - U> <i|), s+ = IT) Ul, 
and s = \[) (t|. Our starting point is an observation 
that if there are M identical atoms, one can straightfor- 
wardly define total spin operator S z = X)j=i s f with 
S ± = ^2j—i sf (j is the index for the atoms), by which 
the atoms represent S = M- — 1, and so on. Keeping 
this in mind, let us consider a coupled array of identical 
cavities in an arbitrary geometry, each of which contains 
M identical single atoms. We employ the Dicke-type 
model, in which every atom in a cavity interacts with 
the cavity mode with the same coupling strength [24j . 
Let us first consider a simple depicted in Fig. [1] 

Let us denote by \tp)j k the state \ip) of the fcth atom in 
the jth cavity. In the rotating frame, the Hamiltonian 
reads 



H=Y, [e iAlt ttiA 



eb 1 „iA 2 t 



2 Af + h.c] 



J2 [(e jAl *5iAf + e^ t g 2 kf)a J + h.c] 
3 

j2 ~(Af + Af ) - + a Ai 



(i) 



U,k) 



where A* v = J2kLi(.\ x ) (y\)jk ( X >V = a , b , e ), a j is tne 
annihilation operator for the jth cavity mode, is the 
corresponding detuning, Qj and j are the corresponding 



Rabi frequencies of the classical fields, gj is the corre- 
sponding atom-cavity coupling rate, and J is the inter- 
cavity hopping rate of photons. Both the transitions are 
coupled to the same cavity mode. The transition between 
I a) and \b) is induced by two-photon Raman transition 
using far-detuned lasers. Here, (j, fc) represent nearest 
neighbor pairs. For now, we ignore the spontaneous de- 
cay rate 7 of the atom. 

Before we proceed, it is instructive to write down our 
parameter regime: 



fll 



Ai A 2 ' 
A h Ai - A 2 > Yy^j » J > , 



AI 



5i_ 
Ai 



Ai 



\lu\ » 2 J. 



(2) 

(3) 
(4) 



The condition @ can be fulfilled with conventionally 
used alkali-metal atoms, such as rubidium and caesium. 
For example, one may choose ground hypcrfme lev- 
els \F = l,m F = -1) and \F = 2,m F = -1) of a 87 Rb 
atom to represent \a) and \b), respectively, and use a + - 
polarizcd light, for which g\ > g 2 . The detuning Aj 
is then comparable to the hyperfine splitting between 
the two levels. Although there are multiple excited lev- 
els, their contributions can be summed up and denoted 
by single parameters in what follows. The other condi- 
tions © and (UJ) can be satisfied simultaneously when 

M gj/Aj 3> J/^J~^fgj. For example, our scheme works 

well in case Aj /1000 - - J, which is allowed 

by strong atom-cavity coupling. 

Our regime is chosen so that the excitation of the atom 
or the cavity photon is suppressed (condition (|3|)), while 
the communication between atoms is mediated by virtual 
cavity photons. The first step is to adiabatically elimi- 
nate the excited state using the conventional method, by 
which the effective Hamiltonian is given by 



^ = -El L ( A r+ A f)«k 



Ai 



E[(MiAf + M2 Af) a J+ h.c] 



(5) 



£|(Af + Af)-£ J(a]a k + aja{), 

3 (3,k) 



where fij = 93 ^ 3 . Now let us introduce spin operators in 
a rotated basis 

lt> = ^(|a> + |6)),|i> = ^(|o)-|6))} (6) 

to represent an s = i spin. Note that these are the 
eigenstates of j(\a) (b\ + \b) (a\). The underlying idea is 
to apply the Raman lasers with Rabi frequency ^ con- 
stantly, introducing a fixed amount of energy splitting 
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\u>\ between the two spin states. The total spin is then 
defined in terms of the operators 



SI 



M 

fc=i 



A I 



4 and ^ = £4- 



(7) 



fe=i 



where s% = ±(|T) (T| - |4) (i\) jk , S+ = (|T) (i\) jk , and 
sjk = (U) (tl)jfe- The total spin is given by 5? = (5f ) 2 + 
+ SrSj). If M is even (odd), the atoms repre- 
sent integral (half-integral) spins up to 4p The atomic 



operators are now written as A 



1 1 



V s" s+ - M_<?^ 

a] t = -K* = f + 5 f • A l x - ^ and A i T = ^ 

Substituting these operators, the Hamiltonian in the ro- 



tating frame reads 



H 



£ 



-e^-)*^S7 La 3 - + /i.e. 



(8) 



where A = M-^- and /uf 2 = /ii ± /i2- Note that in view 
of conditions Q and Q, the effective Rabi frequency 

is much smaller than A, |A ± ui\, and \ui\. This 

allows us to make use of the adiabatic elimination once 
more. We extend the method in Ref. (2f| to keep up to 
the third order terms and take only the subspace with no 
cavity photon. Simple algebra as in Ref. [251 ] yields the 
final Heisenberg spin Hamiltonian 



H=Y J [A{S ] f + B{Sff + CSf] 



£ 



[D(SfS* + SjSl) + ESfS%] , 



(9) 



where A 



~ A 2 - 
- 1 2 
Mi 2 | 
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I A- 



E = 2J| 



Xf-u 2 - ' 2 ' " 2>lA+u 1 \\—u>[ >i ■*"' ""I J 

This Hamiltonian already covers a wide range of para- 
metric regimes for the Heisenberg spin model, although 
individual control of the parameters is limited owing to 
their mutual dependency. Interestingly, the Hamiltonian 
also contains the single-ion anisotropy ( Sf ) 2 , which is of 
essential importance in high-spin cases [261 ] , whereas it is 
merely a meaningless constant in the spin-i case. 

Full control of the individual parameters is allowed by 
bringing in more lasers, shown in FIG. [2] in addition to 
the set up of FIG. [T] The classical fields with Rabi fre- 
quency f^3 and which are applied with an additional 
detuning 5 ~ A, make a similar contribution to the effec- 
tive Hamiltonian as those with Qi and ^2- This can be 
reflected in Hamiltonian (jHJl by adding the same terms 
with A and /x 12 replaced by A — d and /Lt 34 , respectively, 




FIG. 2: Additional lasers to get full control of the spin Hamil- 
tonian. These lasers are applied in addition to the set up of 
Fig.fT] 



where 



^34 



gl^3 I 1 

2 U, 



92^* , 1 
[ A 2 



fi 3 ± fi 4 , 

1 f ). The Stark shift 



aT-r)' and ^ 4 = 
I b) (b\ induced by 



2 va 2 1 A2+5J 
another far-detuned laser field (using a different level and 
polarization) results in adding £\ ^-(e luJt S+ + e~ lu)t S~) 
in Hamiltonian ([5]). Recall that in our previous deriva- 
tion, we have adjusted {0, A, A ± uj} so that they are 
distinct in frequency with the similar frequency spacing 
(condition ([1])), thereby causing each summation in the 
Hamiltonian © to contribute independently to the con- 
stants in the final Hamiltonian ([9]). We adjust ui, A, \±oj, 
A — 5, and A — 5±uj in the same spirit. For the ease of pre- 
sentation, let us take a particular situation where lj > 0, 
A = 3a;, and A — 5 = —6a;, although this is not a neces- 
sary condition. We then obtain the same Hamiltonian ([9]) 

with parameters given by A = -r^Lu^ — M34I )> 

j-y 1/1 I I 91 I ^ 1 I -)- I 91 I \ f~i 

= l' ^12 ~ 16 M12I _ 2 ^34 35 ^34 )' k = 



XvlKI 2 - 

'^3 4 | ), and £ 



Note that C 
2 



16 I A*12 1 70 I ^34 1 )' D — A 2 (32 1^12 

I jo- / ' :i 1 X t(2 |/x£j| +2 l^34| ) 

is determined independently thanks to the term |/i z | z /A, 
while other terms are also determined freely. Hence, this 
parameter set covers any anisotropic or isotropic Heisen- 
berg spin models, with the single-ion anisotropy turned 
on or off. 

The ground state of the spin Hamiltonian could be pre- 
pared by the adiabatic method, as described in Ref. 
Although the Hamiltonian ^ looks similar to a ferro- 
magnetic one (D, E > 0), one can also simulate the anti- 
ferromagnetic spin Hamiltonian, since other parameters 
A, B, and C can be adjusted to have any sign. This can 
be easily seen by noting that if the parameters A, B, and 
C are adjusted to be — 1 times those of a desired antiferro- 
magnetic Hamiltonian, the Hamiltonian is equivalent to 
the antiferromagnetic one up to global factor —1, hence 
with an inverted energy spectrum. Consequently, the 
adiabatic preparation, starting from an antiparallel spin 
configuration ||, ends up with the highest energy state, 
which in fact is the ground state of the corresponding 
antiferromagnetic Hamiltonian. 

Although atomic excitation is heavily suppressed, the 
main source of decoherence in our system is the sponta- 
neous decay of atoms. In relation to the effective spin 
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model, the atomic spontaneous decay results in depolar- 
ization of the spins. This effect can be accounted for 
by considering a conditional Hamiltonian He = H — 
lY^jilA^T + 7s^5 b )' where the effective decay rates 

( \q 2 in i 2 \ 

are approximately given by i A = 7 + [aT+5f ) 

and "/' B — 7 ( J^j- + (2^+(5) 2 ) ' assuming other lasers are 
sufficiently detuned and thus make a negligible contri- 
bution to the decay. In particular, if CljS are chosen in 
such a way that the two contributions are balanced, i.e., 
7^ = "f' B = 7', the depolarization is nearly independent 
of the spin state. In this case, the conditional Hamil- 
tonian is approximately given by He = H — iNM-^, 
where N is the number of cavities. Consequently, the 
state of the system at time t may be written as p(t) = 
e- NM ^' t p Q (t) + (I - e- NM i' t ) P M, where p Q (t) is the de- 
sired quantum state evolved by the spin Hamiltonian and 
Pm is the fully mixed state. This property is useful for 
testing condensed-matter theories, since even under de- 
polarization, the quantum nature retained in the coher- 
ent portion Pq{€) could be observed over a time scale ~ 
1/NMj'. One requirement is that the spin-spin coupling 
rate multiplied by (M/2) 2 should be much larger than the 
global decoherence rate. Reminding that the coup ling 
rate is given by - 2J\p jj \ 2 /X 2 ~ {J/M){Q. j /^/Wj2g j f, 

we require 7 < ^( -M* f . Since A 3 > y/M/2gj 

M j zgj 

from condition ([3]), this requirement can be met for a 
moderate N if J > 7 is satisfied along with our previ- 



ous assumption of strong atom-cavity coupling. Note, 
however, that testing Haldane's conjecture for higher- 
spin chains is more demanding, since the lowest exci- 
tation gap is expected, from its asymptotic behavior, to 
decrease rapidly with increasing M, while the spin cor- 
relation length increases rapidly [161 ]. There are vari- 
ous micro-cavity technologies under active development 
which are expected to fall into our regime of strong atom- 
cavity coupling [2^ . such as superconducting microwave 
cavities [2a|, photonic bandgap microcavities [111, and 
microtoroidal cavities [3(3]. These models would be also 
suited to having a fixed number of atoms in a cavity, by 
virtue of the progress in the micro-fabrication techniques. 
For example, coupling two superconducting qubits with 
a single cavity mode in a well-controlled way has been 
demonstrated recently (28|, which suggests the viability 
of the proposed way of engineering spins [3l| . An impor- 
tant point is that contrary to the system-specific ideas, 
our scheme relies on a general model, which would be 
available in a wide range of current or future systems. 
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